Instabilities of a matter wave in a matter grating 
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We investigate the stability of Bloch waves for a Bose-Einstein condensate moving through a 
periodic lattice created by another condensate modulated by an optical lattice. We show that the 
coupling of phonon-antiphonon modes of the two species give rise to a very rich structure of the 
regimes for dynamical instability, with significant differences with respect to the case of a single 
condensate in an optical lattice. We characterize the relative weight of each condensate in the 
mixing and discuss an analytic limit that accounts for the bare structure of the instability diagrams. 
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Bose-Einstein condensates (BECs) in optical lattices 
have proven to be a versatile tool for investigating funda- 
mental issues of quantum mechanics, solid state physics, 
nonlinear systems and superfluidity [UIS]. Already in ID, 
and in the mean field regime, the dynamics of BECs in a 
periodic potential presents nontrivial features emerging 
from the interplay of nonlinearity - arising from atomic 
interactions - and the discrete translational invariance of 
the lattice. Noticeable features are the formation of soli- 
tonic structures [3 and the presence of several kinds of 
instabilities |i|5lini[Zll3[9lll0l[IIl|l2]. In particular, a 
long investigated issue has been the dynamical instabil- 
ity of nonlinear Bloch waves in the flow of a BEC in a 
ID optical lattice. This phenomenon is characterized by 
the exponential growth of arbitrarily small fluctuations 
that may lead to a partial or complete destruction of the 
condensate [H [HI [10] or to the occurrence of a pulsating 
behavior [TTJ |T2], depending on the interaction regime. 
Recently, also the dynamical stability of a two compo- 
nent BECs has been investigated theoretically [TS", TT . 

All the studies performed so far, however, have been 
restricted to the case where the periodic lattice is not per- 
turbed by the condensate flow, since optical lattices can- 
not support excitations. Different scenarios may occur 
instead when excitations of the lattice interact with that 
of the condensate, like electrons and phonons in a crys- 
tal lattice. This condition can be realized, for example, 
by considering a two component BEC system, where one 
of the two condensates is periodically modulated by an 
external optical potential, playing the role of a periodic 
lattice where the other condensate (transparent to the 
optical potential) flows through. This configuration can 
be designed experimentally by suitably tuning over the 
atomic transitions the electric dipole atom-light interac- 
tion that is used to manipulate the atoms. In particular, 
for alkali atoms in case of large detuning and negligible 
saturation, the intensity of the dipole potential is given 

by m 

(2 + PgFmF , l-PgFmF\ 

where / is the intensity of the laser, /S^x^f are the detun- 
ings from the two fine transitions (I^i e D2) and P is 
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FIG. 1: (Color online) Schematic representation of the sys- 
tem considered: in an interacting two component BEC, one 
of the two condensates (BEC2) is periodically modulated by 
an external optical lattice (dotted line), and acts as periodic 
potential through which the other condensate flows (BECl, 
that is not affected by the optical potential). 



the polarization of the light (P = 0, ±1 for linear and cr^ 
polarization, respectively). For ultra-cold atomic mix- 
tures, since the frequencies of the fine transitions change 
with the species, the potential produced by the same light 
acts in different ways with different atomic species, and 
one can take advantage of this feature by choosing the 
wavelength of the laser so that the potential in Eq. ([T]) 
vanishes for one species but not for the other. 

Actually, a double species BEC with tunable inter- 
species interactions has been recently produced and char- 
acterized Such a system represents indeed a very 
promising tool for exploring novel configurations of ul- 
tracold atomic mixtures in optical lattices [16]. 

In this Letter we study the stability of a one- 
dimensional system composed by two condensates with 
repulsive interactions, under the effect of an optical lat- 
tice to which one of the two species is transparent. In 
this situation it is possible to investigate the effects of 
the mutual interaction in the flow of a non-linear matter 
wave (the condensate that experiences no optical poten- 
tial, hereinafter BECl) through a matter grating (the 
one modulated by the optical lattice, BEC2), as depicted 
in Fig. [1] We find that, as far as dynamical instabil- 
ity is concerned, this system exhibits a very different and 
reacher behavior respect to the usual case of a single BEC 
moving in an optical lattice. In particular, the conditions 
for which BECl becomes dynamical unstable (BEC2 be- 
ing at rest in the lattice) resemble the case of instability 
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for a single attractive condensate [TT], with the addition 
of novel instability regions due to the interspecies inter- 
action. Such instability regions mix up with those for a 
single repulsive condensate when also BEC2 is moving 
with respect to the lattice. We map out the different 
regimes for energetic and dynamical instability, charac- 
terizing the growth rates of the unstable modes and the 
relative weight of each condensate in the mixing. We 
also discuss an analytic limit for vanishing lattice and 
vanishing interspecies interaction, that accounts for the 
bare structure of the instability diagrams in terms of the 
coupling of phono n- ant iphonon modes of the two species. 

The above system can be described by two coupled ID 
Gross-Pitaevskii equations 



dx'^ 



^1 



SCOs2(x)+^2l|^lP+^22|^2p 
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the wavefunctions are normalized to unity over a sin- 
gle period. 

To illustrate a specific case, here we consider the case of 
a ^^Rb -^^K mixture, with an optical lattice of wavelength 
^opt — 790 nm and linear polarization. This choice en- 
sures that the condensate of ^''Rb will experience no op- 
tical potential, playing therefore the role of BECl, while 
the '^^K is subjected to the optical lattice, being BEC2. 
We have also fixed the mean densities to 10 particles per 
lattice site, corresponding to the following values of the 
intraspecies couplings: g\\ c:^ 0.021, ^22 — 0.027. How- 
ever, the effects we are going to present can be obtained 
in a wide range of parameters, with different species or 
different lattice periodicity. 

Equation ^ admits stationary solutions in the form 
of Bloch waves ^j(x) = e^^^^^~^^^ (pjik^x)^ where (j)j{k^x) 
has the same period of the lattice, the Bloch wave vector 
k represents the quasimomentum of the condensate (re- 
stricted to the first Brillouin zone [4 ) , and jUj is the chem- 
ical potential of each species. In the following, for inves- 
tigating the stability of these Bloch states, we consider 
three different configurations: A) the BECl propagates 
through the BEC2 at rest {ki = k^ k2 = 0); B) analogous 
to the former, but with the role of the two condensate 
exchanged {ki = 0, k2 = k); C) both condensate are 
moving with the same quasimomentum {ki = k2 = k). 

As for the case of a single condensate, two kinds of in- 
stability are relevant, namely energetic (or Landau) and 
dynamical instabilities, and can be addressed by consid- 
ering the stability of the system against small perturba- 



tion [4]. For the case of energetic instability the exci- 
tations can be decomposed in terms of Bloch waves of 
quasimomentum ^ as [4] 



1= 



(3) 



where the functions u and v have the same spatial period- 
icity of the lattice. They satisfy the following eigenvalue 
equation 



M{q)X = E{q)X , X 



yd) 



(4) 



obtained by considering the second order contribution to 
the energy functional [4 . The matrix M{q) is given by 



M{q) 



that are written here in dimensionless form, with lengths 
in units of the lattice wavevector ks = ^ir/Xopt^ and with 
the height of the lattice s in units of the recoil en- 
ergy Er = h^k'^/2m2. The coupling constants can 
be expressed in terms of the scattering lengths aij and 
of the mean densities in each site of the lattice Ui as 
gij 2'^^^^0 7r^/n^aijX^^^^^^[l + A(l - Sij)]/k%, with 
A = 1712 /mi (notice that gi2 = 921)- With this choice. 
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(2) 



Ci{q - k) 



(5) 



(6) 



and 



dx'^ 



-\-ik] + (j — l)s cos^(x) 



(7) 



(8) 



The energetic instability occurs when the condensate fiow 
corresponding to a Bloch wave is no more a local min- 
imum of the energy functional, that is when the Mk{q) 
matrix has negative eigenvalues [4 . 

For the case of dynamical instability we look for time- 
dependent fiuctuation of the form 

^ /Ti a \ / 



yielding the following Bogoliubov equations [H [131 E] 

I]3M(q)X = ^uqX , E3 = 'J ) (9) 

where (73 is the third Pauli matrix. Dynamical in- 
stability occurs when some of the eigenfrequencies uokq 
get a nonzero imaginary part, indicating an exponential 
growth of the corresponding modes, that rapidly drives 
the system away from the steady state in the dynamical 
evolution. 

The typical stability diagrams in the k — q plane are 
shown in Fig. [2] (for the case A and B) and in Fig. [s] 
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FIG. 2: (Color online) Configurations A (left) and B (right): stability diagrams as a function of the quasimomenta k and q of 
the condensates and the excitations, for different values of s and gi2. Regions in gray (color) scale indicate where the system is 
dynamically unstable. The gray (color) scale corresponds to the growth rate of the unstable modes {(xlm{ookq))- Shaded area: 
energetic instability. Solid lines: coupling between phonon-antiphonon modes in the limit s ^ 0, gi2 (see text). 



(case (7), for different values of the lattice intensity s 
and interspecies interaction gi2. In all three cases, the 
energetic stability behavior is similar to that of the single 
species [4], whereas the regions of dynamical instability 
show a very different and richer behavior, owing to the 
coupling between the two condensates, as we will explain 
below. Moreover, the fact that for large enough repulsion 
between the two species the system becomes unstable 
also for low k values (see the bottom row diagrams) is 
related to occurrence of a phase separation instability 
[la [H [IHl mEO] for > 7^11^22 (here the critical 
interaction strength is — 0.018). 

Let us now analyze in more details the structure of 
the dynamical unstable regions. Since these instabili- 
ties are determined by the mixing of both components, 
it is instructive to know to which extent they con- 
tribute. To measure this effect here we define Vi(k^q) = 

I'^i* ^ introduce the following 

quantity 



r]{k,q) 



Vi{k,q)-V2{k,q) 
Vi{k,q)^V2{k,q) 



(10) 
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FIG. 3: (Color online) Similar notations of Fig. [2]but for con- 
figuration C. Here the gray (color) scale corresponds to the 
values of r](k, q) in the dynamically unstable regions, measur- 
ing the relative weight of each species. 



that is positive or negative when the major component 
corresponds to the BECl or BEC2 respectively, and close 
to zero when the contribution of the two BECs is almost 
the same. This indicator is shown explicitly for the con- 
figuration C in Fig. |3) since this is the case that presents 
the most complex mixing behavior (due to the fact that 
both condensate quasimomenta are non vanishing). In 
particular, it emerges that the instability at low values 
of k gets the primary contribution from BECl, whereas 
that at higher k is dominated by BEC2. Instead, for the 
configurations A and B of Fig. [2]7^ is mostly positive (A) 
or negative {B). 



Further insight on the origin of the different instabil- 
ity regions can be gained by considering the limit 5^0 
and gi2 0. In this case the full matrix TisM{q) is block 
diagonal, the two species are decoupled, and the eigenval- 
ues - that will be indicated as cc;!^ {q) (the ± refers to the 
positive and negative sectors of the spectrum, i = 1,2) 
- have a simple real analytic expression in all the three 
configurations [4^. In this limit therefore the system is 
dynamically stable. However, as soon as s 7^ and 
9i2 7^ 0, a resonant coupling between phonon and anti- 
phonon modes may cause the appearance of an imaginary 
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FIG. 4: (Color online) Analytical solutions for the coupling 
of the phonon and anti-phonon modes in the limit s ^ 0, 
9\i for the three configurations A,B and C. 

component in the spectrum, and therefore dynamical in- 
stability [21 . Notably, in all the three cases, there are 
several ways to couple phonon and antiphonon modes, 
due to the fact that besides the usual intraspecies cou- 
pling [4 here we can also have interspecies resonances, 
namely — 2) = u;'^_{q)^ ^+(^ ~ 2) = u;^_{q)^ and 

uj'^±{q) = The curves corresponding to the possi- 

ble resonances lying in the range k^q G [0,1] are reported 
in Fig. [4] for each of the three configurations A,B or C. 



They are also reported in the top left panels in Figs. [2]- 
[3] Then, as we move away from the analytic limit, both 
by increasing the interspecies interactions gi2 or the lat- 
tice intensity 5, the strong mixing between the different 
modes eventually causes a broadening and a deforma- 
tion of the instability regions. Remarkably the seeding 
of most of the instabilities is the coupling of a phonon 
(antiphonon) mode of the species 1 with an antiphonon 
(phonon) mode of the species 2. This situation is com- 
pletely different with respect to the case of a single con- 
densate in an optical lattice, since the latter cannot sup- 
port excitations. 

Summarizing, we have investigated the stability of non- 
linear Bloch waves in a one-dimensional two component 
EEC system, in a configuration where one of the two 
condensate is modulated by an external optical lattice, 
effectively realizing a matter wave periodic potential for 
the other condensate. Such a system presents a very rich 
structure of the regimes of dynamical instability, origi- 
nating from the resonant coupling of excitations of the 
two BECs. Different scenarios can be obtained by tun- 
ing the interspecies interaction and the relative velocity 
between the two condensates and/or the external optical 
lattice, allowing for different ways to match the resonance 
of phonon antiphonon branches of the two species. All 
these features make multicomponent ultracold atoms in 
optical lattices a very promising system for studying the 
interplay of nonlinear ity and periodicity in novel regimes. 
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Owing to the properties of the matrix T>3M(q), complex 
eigenvalues can appear only in pairs, so that the onset of 
instability is signaled by the appearance of a pair of real 
degenerate eigenvalues [4^- 



